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jy^ • We explain the existence of excess events near the endpoints of the double 

beta decay (pp) spectra of several elements, using the neutrinoless emission of 
massless Goldstone bosons. Models with scalars carrying lepton number —2 
are proposed for this purpose so that ordinary neutrinoless 00 is forbidden, 
and we can raise the scale of global symmetry breaking above the 10 keV scale 
! needed for observable emission of conventional Majorons in 00. The electron 

spectrum has a different shape, and the rate depends on different nuclear 
matrix elements, than for the emission of ordinary Majorons. 
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Double beta decay {(3(3) is an extremely rare process in which two neutrons in a 
nucleus simultaneously decay into protons. It has now been observed in seven elements: 
76 Ge, 82 Se, 100 Mo, 128 Te, 130 Te, 150 Nd and 238 U. One reason (3/3 has received so much 
attention is the possibility of detecting small Majorana masses for the neutrinos if the 
neutrinoless mode /3/3oi/ should be observed. The amplitude for (3(3^ is proportional to an 
effective Majorana mass, m e fj = ^2 @i m ii where 9i is the mixing angle between v e and the 
ith mass eigenstate. Experimentally it is known that 



X '■ rn eff < 1 eV. (1) 

sh ; 

It is also possible that f3(3o u proceeds by the annihiliation of the virtual neutrinos, 1 ' 2 
through a vertex 

gifUe^Ue, (2) 

where a scalar particle f is emitted. This can naturally occur if lepton number is broken 
spontaneously, so that neutrinos couple to a massless Goldstone boson called the Majoron. 
The sum-energy electron spectrum for this process, which we refer to as (3 (3m, is skewed 
toward higher electron energies than that for fifov, because it is a three-body decay. 
Previous experiments searching for this mode obtained a limit 3 of g < 2 x 10 -4 , following 
indications of a positive signal at the 3 x 10 -4 level. 

It is therefore intriguing that the UC Irvine group has for the last few years been seeing 
excess events near the endpoint in three different elements: 82 Se, 100 Mo and 150 Nd. 5 A 
similar excess is also plainly visible in the 76 Ge spectrum of Avignone et al. 6 In all four 
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elements, the excess events become evident starting from 0.5 to 1 MeV below the endpoint, 
and constitute between 2 and 3 % of the total signal. Because the expected signal from 
P^2v is negligible near the endpoint, the chance of these events being due to statistical 
fluctuations is only 1 in 10 5 . 

Let us write the partial rate for (3 /3m as 

dT 0vM = 2n- 5 (g A G F ) 4 g 2 eS \M\ 2 dn, (3) 

where g e ^ is a model-dependent effective coupling, Ai is the usual combination of Gamow- 
Teller and Fermi nuclear matrix elements, and <ifi is the phase space. We have analyzed 
the data by matching the number of excess events above some threshold energy where 
visually the anomalies seem to begin for each element. For 238 U and the Te isotopes, only 
the total rate is observed, so we omit thresholds for these elements to see how large a 
Majoron coupling would be needed for (3 (3 m to saturate the total observed signal. 

The results are shown in Table 1. We find that the required values of g e g are in the 
range 5 x 10 -5 to 2 x 10 -4 except for 100 Mo, which needs a somewhat larger coupling. 
This presents a problem for the ratio of the rates of 130 Te to 128 Te decays, which has 
recently been measured 7 to be (2.41 ±0.06) x 10 3 . The (3 (3m prediction implied by Table 1 
(and the phase space calculations not shown there) would be much smaller, 93, assuming 
<7eff = 1 x 10 -4 so that the endpoint anomalies could be explained. 



Table 1: couplings for emission of ordinary Majorons in double beta decay, assuming total rate of 
geochemically observed decays was saturated by (3(3 M , and using the nuclear matrix elements of 
ref. [8]. 



Element 


76 Ge 


82 Se 


100 Mo 


150 Nd 


128 Te 


130 Te 


238-y 


9eS 


lxlO" 4 


8xl0" 5 


4xl0" 4 


8xl0" 5 


4xl0" 5 


lxlO" 4 


2xl0" 4 



Even if we ignore the Te problem, there is a more serious difficulty: no Majoron 
models exist which seem able to give such a large coupling. The favored theory for (3(3m 
was the triplet Majoron model, 2 ' 9 but it has been ruled out by LEP's measurement of 
the invisible Z width. Alternatively, the singlet Majoron model 10 has couplings which are 
highly suppressed, since generically 

9eS ~ rnes/v, (4) 

where v is the scale of lepton number breaking. The most natural choice for v is of order 
the weak scale, leading to g e Q ~ 10 -11 . But we need v < 10 keV to explain the size of the 
(3(3 anomalies, given the limit (1) on m v . This unnaturally small value requires extreme 
fine tuning to keep it so far below the weak scale.* However, suppose we could remove 
the relation between v and m e ff by maintaining lepton number as an exact symmetry, 
thus insuring m e ff = to all orders. Perhaps we would be able to have a higher scale of 
breaking for whatever symmetry it is that gives us the Goldstone bosons. But if lepton 
number is conserved, then the Majorons must carry —2 units of L. We will call these 
charged Majorons since they carry the global lepton charge. 



* Moreover, even if v ~ 10 keV, one must introduce some neutrinos heavier than the 100 MeV 
scale of the nucleon Fermi momentum; otherwise the mixing of the light neutrinos is ineffective 
and g e {t becomes the bare coupling of the Majoron to v e , which is zero for singlet Majorons. 



This situation reminds us of the standard model, where a would-be Goldstone boson, 
the longitudinal component of W~ , has electric charge. It is therefore quite easy to con- 
struct a model of charged Majorons; we just mimic the standard model by introducing an 
SU(2) s xU(l)i/ global symmetry that acts on the leptons and on some new sterile neutri- 
nos and Higgs multiplets. Although we could give a renormalizable model of this sort,* it 
is clearest to present an effective Lagrangian with the dimension-five operators 

f (L e tf)(l + 75)0^) + £(* f *)(l + 75)(0 f s) + h.c, (5) 

where <p = T2<p* ■ These particles transform as L e ~ (1,1,2,-1/2), s ~ (2,0,1,0), <p ~ 
(2,1,1,0) and H ^ (1,1,2,-1/2) under the SU(2) s xU(l) L >xSU(2) L xU(l) y symmetries. 
When (j) gets a VEV v, a Dirac mass term between s± and S2 results. v e is mostly a massless 
state v, but mixes with a heavy Dirac neutrino ^ with mixing angle 9 = tan -1 (g (H) /hv). 
One finds that the charged Majorons couple to the neutrinos via 

I 

— (sin Qv - cos dW) (1 - 7 5 )^* c + h.c. (6) 
Av 

They have charge —2 under the residual lepton symmetry which is given by L = L' plus 
the diagonal generator of SU(2) S . 

We note that the sum-energy electron spectrum for emission of charged Majorons 
differs from that of ordinary Majorons. Consider singlet Majorons: in the language where 
they are derivatively coupled (v ~ 1 u0ifj^u), there exist extra Feynman diagrams with Ma- 
joron Brehmsstrahlung off the electron lines via the vector coupling v~ l e$(pe, while the 
neutrinos annihilate through their Majoron mass. These survive at zero Majoron momen- 
tum q because the internal electrons are going on shell as q — > 0. But these graphs don't 
exist for charged Majorons because there are no Majorona masses; hence the amplitude 
vanishes as q — > 0, and consequently the sum energy spectrum of the electrons has a shape 
intermediate between that of fifav and generic (3(3m spectra. Our spectrum behaves like 
(Q — E) 3 near the endpoint energy Q, compared with (Q — E) 5 for P@2v and (Q — E) for 
ordinary (3 (3m- The different spectra are shown in Fig. 1. 

Furthermore (3(3m with charged Majorons depends on different nuclear matrix ele- 
ments than the usual ones, because the amplitude contains an extra factor of the neutrino 
momentum p. The leptonic part has the Lorentz structure $,4\lv eC L , an d for kine- 
matic reasons only the spatial components of p contribute significantly. Since they have 
odd parity, but the nuclear transition is + — > + , we require nuclear operators with odd 
parity also, unlike the Gamow- Teller and Fermi terms. Such operators come from the 
recoil corrections to the hadronic weak currents, as well as p-wave Coulomb corrections to 
the electron wave functions. For example, we have matrix elements of the form 

M R = (9v/9a) (p 2 )' 1 (r ■ (D n x <r m + a n x D m ) JU(r)^ , (7) 

using the notation of Doi et al. 11 or Tomoda. 12 



* For example, introduce sterile neutrinos Ni and A^2 having U(l) L i charges +1 and —1 re- 
spectively, and include all their allowed couplings to L e , H, (f> and s. Integrating out Ni results in 
the effective theory. Even if Ni is no heavier than the remaining particles, the resulting prediction 
is quantitatively similar to that of the renormalizable theory. 



Fig. 1: Two-neutrino and Majoron emitting electron sum energy spectra for Ge. 



We find our process is roughly comparable to the ordinary Majoron process with a 
coupling of 

/,, /0 ^2 (qp)M 2 M new 
g«~{M/2v)0 {{p2)+M2)2 M , (8) 

where M is the mass of the heavy sterile neutrino, 9 is its mixing angle to v e , and (qp) is of 
order 100 MeV 2 since the Majoron only has energy q ~ 1 MeV, whereas the virtual neutrino 
momentum is ~ 100 MeV. -M new denotes the new nuclear matrix elements associated with 
charged Majorons. Note that M ~ v is also the scale at which our global symmetry breaks. 
To get <7eff as large as possible we need 9 ~ 0.1, which forces us to take M > 500 MeV due 
to the constraints from peak searches in meson decays. 13 ' 14 With these values, we would 
also need the new nuclear matrix elements (some of which have not been computed, to 
our knowledge) to be 25 times larger than the usual ones. Another possibility is that the 
coupling producing the heavy neutrino mass is very strong (M > v), and that it decays 
so quickly into light neutrinos and Majorons that it is really a rather broad resonance. 
Then its mass might be closer to 100 MeV and yet escape notice in it and K decay peak 
searches, even with a large mixing angle. 

More exactly, the rate for charged Majoron emission is, 

dT M = 2n- 5 (g A G F 9) 4 (M/2v) 2 \M new \ 2 dn new . (9) 

The detailed forms of -M new and <ifi new are given in ref. [15]. Some of the contributions to 
.M n ew, such as the one displayed in eq. (7), have been calculated. In Table 2 we compare 
the value of -M ne w needed to account for the anomalous events with the computed value 16 
of M. R , one of the actual contributions to A4 new . (For this comparison we have assumed a 
large mixing angle and strong coupling, as indicated in the table caption.) The agreement 
is remarkably good except for 128 Te. However, the absolute rate 128 Te is known less well 
than the ratio to 130 Te, and our model predicts r( 130 Te)/r( 128 Te)= 770, in much better 
agreement with the experimental value quoted above than in the singlet Majoron model. 

In summary, we have presented a model of Goldstone boson emission in double beta 
decay which may be able to explain excess events near the endpoints of several elements. 
The boson has lepton number L = — 2 and the light neutrinos remain exactly massless so 



Table 2: matrix elements for charged Majoron emission in double beta decay. We give the values of 
the matrix elements which are needed to explain the data (assuming all events were f3/3 M for 238 U 
and Te) for 9 2 (M/2v)M 2 /((p 2 )+M 2 )=3xlO~ 2 , and values of some representative matrix elements 
that have been calculated. 16 



Element 


76 Ge 


82 Se 


100 Mo 


150 Nd 


128 Te 


130 Te 


238tj 


TV! n ew needed 


1.1 


0.99 


1.4 


1.3 


0.32 


0.48 


0.61 


\M R \ (Muto et al) 


1.1 


0.95 


1.1 


1.3 


0.92 


0.78 


? 



that ordinary neutrinoless double beta decay is forbidden. The electron spectrum for these 
double beta decays with "charged Majoron" emission is less peaked toward high energies 
than that for ordinary Majorons. The model requires a sterile Dirac neutrino with a large 
(~ 0.1) mixing angle and hence (due to experimental constraints) a mass M ~ 500 MeV. 
Moreover, the (3(3m rate in this case depends on unknown nuclear matrix elements A4 new , 
different from those appearing in the usual amplitudes for 1313m, which must turn out to 
be larger than the usual ones in order for the rate to be as large as is seemingly observed. 
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